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Abstract. We present the 2 — 200 /im Infrared Space Observatory (ISO) spectrum and 3.8 — 20 fim ISAAC and 
TIMMI2 images of the extreme OH/IR star IRAS 16342-3814. Amorphous silicate absorption features are seen 
I . at 10 and 20 /im, together with crystalline silicate absorption features up to almost 45 fim. No other OH/IR 

star is known to have crystalline silicate features in absorption up to these wavelengths. This suggests that IRAS 
■^j- ' 16342-3814 must have, or recently had, an extremely high mass-loss rate. Indeed, preliminary radiative transfer 

^*) , calculations suggest that the mass-loss rate may be as large as 10 -3 Moyr -1 . The 3.8 /xm ISAAC image shows a 

t-H ■ bipolar reflection nebula with a dark equatorial waist or torus, similar to that seen in optical images taken with the 

Hubble Space Telescope (HST). The position angle of the nebula decreases significantly with increasing wavelength, 
suggesting that the dominant source of emission changes from scattering to thermal emission. Still, even up to 
20 fj,m the nebula is oriented approximately along the major axis of the nebula seen in the HST and ISAAC 
images, suggesting that the torus must be very cold, in agreement with the very red ISO spectrum. The 20 /im 
image shows a roughly spherically symmetric extended halo, approximately 6" in diameter, which is probably 
due to a previous phase of mass- loss on the AGB, suggesting a transition from a (more) spherically symmetric 
to a (more) axial symmetric form of mass-loss at the end of the AGB. Using a simple model, we estimate the 
maximum dust particle sizes in the torus and in the reflection nebula to be 1.3 and 0.09/im respectively. The size 
of the particles in the torus is large compared to typical ISM values, but in agreement with high mass-loss rate 
objects like AFGL 4106 and HD161796. We discuss the possible reason for the difference in particle size between 
JL " ■ the torus and the reflection nebula. 
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1. Introduction 10~ 5 M Q yr~ 1 . The SW rapidly strips the star of virtually 

its entire envelope. The high mass-loss during the AGB 

As a Main Sequence star with an initial mass between phage (thg ^ AQB phage , agtg gQme 1q6 yp) createg & 

1 and about 8 M ages, it eventually evolves into an circumstellar enve lope which may completely obscure the 

Asymptotic Giant Branch (AGB) star. An AGB star ^ when the gw ph&ge gtopgj ^ centr&1 gt&r gtartg tQ 

is a luminous cool giant star in which a dust driven , . , . , , , , . ., ... 

° snrmk m radius at constant luminosity, getting hotter over 

wind expels matter from its envelope at a high rate ^ Thfi gt&r ig nQw gaid tQ be ^ the post _ AGB phase 

if the* AGB t^^J^^t^^J^'^Z The ° ld circumstellar envelope, which initially obscures 

as the GB star evolves. This eventually ends in an ^ e central star, becomes optically thin at optical wave- 

episode of extremely high mass-loss, the so-called super- lengthg after _ 1QQ _ 1QQQ yr bec&uge - t ig gtm expanding 

wind (SW) phase, in wh ich the mass-loss rate, M, exceeds Depending on the rate at which the central star evolyeS) 

Send offprint requests to: C. Dijkstra (dijkstra@astro.uva.nl) the star ^ start to photoionize the material that was 

* Based on observations with ISO, an ESA project with in- e J ected dunn S the AGB P hase ' creating a planetary neb- 

struments funded by ESA Member States (especially the PI ula (PN). In this case the phase between the end of the 

countries: France, Germany, the Netherlands and the United AGB and the birth of the PN is called the proto-planetary 

Kingdom) with the participation of ISAS and NASA. The SWS nebula (proto-PN) phase, 
is a joint project of SRON and MPE. Also based on observa- 
tions obtained at the European Southern Observatory (ESO). 
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In this paper we study the proto-PN IRAS 16342-3814 
(hereafter: IRAS 16342; also OH 344.1+5.8 or "the water- 
fountain nebula"). The evolutionary status of IRAS 16 342 
as a proto-PN is examined in detail by ISahai et alJ 
l)l999t hereafter: SAH99). They report the presence of 
an asymmetrical bipolar reflection nebula in Hubble Space 
Telescope (HST) images at 0.55 and 0.80 /im. The two 
lobes of the nebula, where the eastern lobe is significantly 
redder than the western lobe, are separated by a dark 
equatorial waist. The nebula is illuminated by starlight 
escaping through the polar holes in this 'optically thick, 
dense, flattened cocoon of dust'. The central star is com- 
pletely obscured by the dust cocoon (hereafter called 'the 
torus'). The individual symmetry axes of the lobes do not 
coincide, although their orientation and location with re- 
spect to the dark waist are centrally point symmetric. OH 
maser line data at 1612, 1665 and 1667 MHz overlaid on 
the 0.80 fim image suggest that the nebular axis is tilted 
such that the western lobe is nearest to us. Faint halos 
are seen around the lobes. The halo regions are illumi- 
nated by scattered light from the lobes, and possibly by 
heavily reddened, direct starlight (SAH99). 

IRAS 16342 has been ob served as part of a mid- 
infrared llMeixner et alJll999|) and optical imaging survey 
l(Ueta et all200(t . In the mid-infrared survey IRAS 16342 
was observed at 9.8, 12.5, 17.8 and 20.6 fim. At these 
wavelengths the star was classified as a core/elliptical 
source, i.e. a source with an unresolved, very bright core 
that is surrounded by a lower surface brightness elliptical 
nebula. At all mid-infrared wavelengths the source was 
found to extend northeast-southeast with the same posi- 
tion angle for all wavelengths. 

Early studies on IRAS 16342 revealed the pres- 
ence of extremely high velocity outflows in both H2O 
(~130 km s -1 ) and OH (^70 km s" 1 ) maser lin e emission 
i Zuckerman Lolll 987t lbikkel fc Morrislll 988h . Likkel & 
Morris interpreted the H2O masers as streams or clumps 
of molecular gas forced out at the polar axes, and believed 
the OH masers to origin ate from lower velo city material 
at intermediate latitudes. iLikkel et al 1 lll992t suggested a 
binary star model for IRAS 16342, as this could account 
for both the bipolarity and the high outflow velocities of 
the masers. 

Radiative transfer ca lculations on IRAS 16342 
have been perform ed by iGuertler et al ] lll996ft and 
IZiilstra et al I l|200lh . The latter authors found an inner 
radius of 1" and an outer radius of 1.5" for the torus, 
in agreement with the HST observations. The mass of 
the torus was found to be 0.1 M Q . The radiative transfer 
calculations by both groups suggest the presence of a 
(very cool) outer shell. 

In this paper we present the 2 — 200 /im Infrared Space 
Observatory (ISO) spectrum and infrared ISAAC and 
TIMMI2 images (3.8 to 20 pm) of IRAS 16342. The aim 
of this paper is to study the mineralogy, morphology and 
mass- loss history of IRAS 16342. The paper should also 
add to our understanding of proto-PNe in general and 
to the injection of matter into the ISM by AGB stars. 



Table 1. The IRAS fluxes for IRAS 16342 compared to 
those derived from the ISO data. From left to right we list 
the wavelength; the Point Source Catalog (PSC) IRAS 
flux; the IRAS flux derived from the ISO spectrum. For 
more details see Sect. 12.11 



A 


PSC IRAS 


ISO flux 


H 


flux (Jy) 


(Jy) 


12 


16.2 


12.0 


25 


199.8 


170.3 


60 


290.2 


292.5 


100 


139.4 


165.3 



The latter is important since AGB stars are among the 
main 'factories' of interstellar dust in galaxies. In Sect. 2 
the data and data reduction procedures are described. In 
Sect. 3 we present and discuss the ISO spectrum. Based 
on the spectrum, we study the mineralogy and mass-loss 
history of IRAS 16342. The ISAAC and TIMMI2 images 
of IRAS16342 are presented in Sect. 4 and will be used to 
study the morphology of IRAS 16342 and set limits on the 
grain sizes in the torus and the lobes. The conclusions are 
listed in Sect. 5. 

2. Data and data reduction 

2.1. ISO data 

The 2 — 200 /im spectrum was obtained with ISO (Kessler 
et al. 1996), using the Short Wavelength Spectrometer 
(SWS; Ide Graauw et al.lll996|) for the 2 - 4 5 region 
and t he Long Wavelength Spectrometer (LWS: lciegg et all 
1996) for the 45 — 200 /im region. The two spectra were 
reduced independently. For both the SWS and LWS re- 
duction, we used pipeline version OLP 10 as the starting 
point. In both cases, individual subbands were cleaned 
from glitches, flatfielded (to the mean of all detectors), 
sigmaclipped (using the default values a = 3 for SWS 
and a = 2.5 for LWS) and rebinned (A/AA = 500 for 
SWS and AA = 0.068 /im for LWS). The subbands were 
scaled (to band 1 and/or 2 for SWS and detector 5 for 
LWS) to obtain a continuous spectrum. Finally, the LWS 
spectrum was scaled to the SWS spectrum to get one con- 
tinuous spectrum. Without scaling, the flux levels of the 
two spectrographs matched at 45 /*m to within 12%. For 
LWS there was an off source spectrum available to cor- 
rect for interstellar flux contributions. The flux levels in 
this spectrum were low (<10 Jy), showing that the on- 
source spectrum has no significant background contribu- 
tion. We have compared our reduction with the available 
IRAS data. Between 14 and 25 /im, the ISO spectrum 
overlaps the IRAS Low Resolution Spectrum (LRS) to 
within 5%, while between 12 and 14 /im the IRAS LRS 
flux is about 25% higher than the ISO flux. The IRAS 
Point Source Catalog (PSC) fluxes match those derived 
from the ISO data to within 26% (see Table [TJ. 



C. Dijkstra et al.: The mineralogy, geometry and mass-loss history of IRAS 16342-3814 



3 



2.2. ISAAC data 

The IRAS 16342 3.8 fim (L-band) image was obtained 
with the Infrared Spectrometer And Array Camera 
(ISAAC) 1 at the Very Large Telescope (VLT) at Paranal, 
Chile, on August 6, 2000, using the Santa Barbara 
Research Center Aladdin 1024 x 1024 pixel InSb array. The 
pixel scale of this detector is 0.071 " per pixel. The total 
field of view is 72" x 72". The L-band filter has a central 
wavelength of 3.78 /jm and a width of 0.58 fim. The op- 
tical seeing during the observations varied between 1.2" 
and 2.0". 

The image was obtained using a 10" chop throw (per- 
pendicular to the long axis of the source), while nod- 
ding was performed. In addition to the chopping, the 
source position was also jittered. The total integration 
time was 5.0 minutes, not considering the negative im- 
age in the chopped image. Only the positive image was 
used, since the PSF of the negative image was systemat- 
ically wider. We used two photometric calibration stars: 
Y3501 and HR6736. The L-band magnit udes of these stars 
wer e measured by MSSSO ph otometry l|McGregorlll994|) 
and Ivan der Bliek et al] l|l996t) , respectively. As a refer- 
ence for the point spread function (PSF), a nearby star, 
USNO 450-24728679 (m L = 7.0 magn) was used. The 
obtained PSF was used to deconvolve the ISAAC image 
of IRAS 16342. The FWHM of the PSF was calculated to 
be 0.35". The reduction was done using eclipse package 
version 4.1.2 jPevillardl Il997|) . and DAOPHOTO under 
IDL. 

2.3. TIMMI2 data 

IRAS 16342 images at 7.9 (narrow band), 9.8 (narrow 
band), 10.6 (broad band) and 20.0 fim (broad band) 
were observed with the Thermal Infrared Multi Mode 
Instrument (TIMMI2) at the ESO 3.60m telescope at 
La Silla, Chile, on August 11, 2001, using a Raytheon 
240 x 320 pixel AsSi array. The pixel scale of this detector 
is 0.202" per pixel. The recorded optical seeing during the 
observations varied between 0.92" and 2.0". 

The images were obtained using a 15" chop throw in 
the north-south direction while nodding was performed. 
As a reference for the PSF, HD 124897 (= a Boo) and 
HD 178345 were used. We combined individual chopped 
images into one image to improve the signal-to-noise. This 
was done by shifting the positive and negative images in 
each chopped frame to a prescribed position and then 
adding or subtracting them. Frames of poor quality were 
omitted. Finally, we deconvolved the resulting image, us- 
ing either HD 124897 or HD 178345 as PSF. Due to non- 
photometric weather during the observations, the absolute 
flux calibration of the TIMMI2 images proved difficult and 
the results uncertain. Therefore, further attempts to do 
an absolute flux calibration on the images were not made. 
Moreover, we are mainly interested in the morphology of 

1 http:/ /www. eso.org/instruments/isaac/index. html 



IRAS 16342 at different wavelengths, which does not re- 
quire an accurate (absolute) calibration of the flux. We 
will use relative fluxes in the individual images only. 

3. The ISO spectrum 

3.1. Description of the spectrum 

In the top panel of Fig. ^ we present the 2 — 200 fim ISO 
spectrum of IRAS 16342 (solid line). A possible contin- 
uum for the spectrum, which we will discuss below, is 
also given (smooth solid line). The second panel shows 
the spec trum of the OH/IR st ar AFGL 5379, studied in 
detail bv lSvlvester et all (|l999h . The continuum for AFGL 
5379 is taken from that paper. The 10 and 20 /im amor- 
phous silicate features are clearly visible. These absorption 
features are observed in many AGB stars, their strength 
depending on the optical depth of the dust, which is di- 
rectly related to the mass-loss rate. The absorption fea- 
tures become stronger, and the spectrum redder, as M 
increases (e.g. Bedijn, 1987). IRAS 16342 also shows the 
amorphous silicate features. The one at 20 /im was al- 
rea dy detected in the IRAS LRS spectrum and fitted 
by IZiilstra et ail l|200lh . Compared to AFGL 5379, the 
amorphous silicate features in the spectrum of IRAS 16342 
might seem to be absent at first glance! The spectrum of 
IRAS 16342-3814 is that of an extremely reddened OH/IR 
star. The third panel in Fig. ^ shows the spectrum of 
HD 161796. HD 161796 is ano ther post- AGB object w hich 
has been studied in detail bv lHoogzaad et alJ l|2002|) and 
iMeixner et"aT] ll2002t). The contin uum for HD 161796 is 
taken from Hoogzaad et alJ l)2002j) . 

A first inspection of the spectrum of IRAS 16342 sug- 
gests the presence of the characteristic features of crys- 
talline silicates and crystalline water ice, though up to 
45 fim they appear in absorption instead of in emission. 
To illustrate this, the continuum shown in the top panel 
of Fig. [T] has been drawn accord ingly (smooth solid line). 
Following iMolster et all <|2002bj) . we represented the con- 
tinuum by a spline-fit and required it to be smooth and 
to maximize the continuum flux in both F u and F\. We 
emphasize that this continuum is used only to enhance 
the visibility of sharp features in the spectrum and has no 
physical meaning. In the lower panel of Fig. ^ we plot 
the continuum subtracted spectra of IRAS 16342 (solid 
line) and HD 161796 (dashed line). The continuum sub- 
tracted spectrum of HD 161796 is dominated by crys- 
talline silicate and cry stalline water ice emission features 
(Ho ogzaad et al. 2002). A comparison of the suggested 
absorption features in IRAS 16342 with the emission fea- 
tures in HD 161796 shows an impressive agreement, indeed 
identifying the IRAS 16342 features as due to crystalline 
silicates and water ice, and justifying our original choice 
of the continuum. At wavelengths longwards of approxi- 
mately 45 fim the features in the spectrum of IRAS 16342 
are again in emission. Close to the peak of the SED, where 
the transition between absorption and emission occurs, the 
shape of the spectral features is less certain. 
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Fig. 1. The 2 - 200 /mi ISO spectra of IRAS 16342-3814 (top panel), the OH/IR star AFGL 5379 (second panel) and 
the post-AGB star HD 161796 (third panel). Adopted continua for the spectra are given (smooth solid lines). The 
lower panel shows the continuum subtracted spectra of IRAS 16342 (solid line) and HD 161796 (dashed line). For a 
discussion on sec Sect. 13. ll 




Wavelength [/xm] 



3.2. The mineralogy of IRAS 16342 

Tabic lists the wavelength positions and identifications 
of the features found in the spectrum of IRAS 16342. 
We identify amorphous silicates, forsterite (Mg 2 Si04), 
diopside (MgCaSi20e), crystalline water ice, and possibly 
clino-enstatite (MgSiOa) as important solid state compo- 
nents of the circumstellar matter around this object. We 
also find [01] and [CII] lines, indicating low density gas. 
The off source spectrum suggests that the [CII] line is 
largely interstellar, although a circumstellar contribution 
(up to 25% in peak strength) may still be present. The 
[01] line seems to be interstellar. After subtraction of the 
off source spectrum, in which the [01] line is clearly seen, 
there is no clear evidence for a residual circumstellar com- 
ponent. 

Fig. shows part of the continuum subtracted spec- 
trum of IRAS 16342 in more detail, together with the 



mass absorpti on coefficients for forsterite, clino-enstatite 
and diopside l|Koike et alJ 1999). Also shown is the ab- 
sorbance of cryst alline water ice measured at 50 K 
l|Smith et alJll994h . The wavelength positions of the iden- 
tified forsterite, diopside and crystalline water ice features 
are indicated with short vertical lines. We also indicate the 
positions of clino-enstatite bands. Due to blending with 
forsterite, and given that the bands are fairly weak, we 
cannot be certain that clino-enstatite is present. However, 
in most evolve d stars both forsteri te and enstatite are 
found together l|Molster et alJl2002aj) . and we suggest this 
is also the case in IRAS 16342. The 18.5 fxm amorphous 
silicate feat ure, the [OI] line and the unidentified 48 //m 
feature (e.g. iMolster et alJEo02a|) are also indicated. 

The solid state materials found in the circumstel- 
lar environment of IRAS 16342 are also found around 
other (high mass- loss), oxygen rich, evolved stars (e.g. 
iKoike et alJ l2000 HMolster etaTll2002ah . The very red ISO 
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Fig. 2. The 15 — 55 /mi and 40 — 100 /xm continuum sub- 
tracted spectrum of IRAS 16342 together with the ma- 
terials identified in it. The mass absorption coefficients 
for forsterite (F), clino-enstatite (E) and diopside (D), to- 
gether with the absorbance spectrum of crystalline water 
ice (H20) are also shown. The continuum in the mass ab- 
sorption coefficient and absorbance curves has been sub- 
tracted in order to enhance the spectral features of the ma- 
terials. Identifications of spectral features are given where 
possible. Note that arbitrary scaling and offset values have 
been applied to the laboratory data in order to reveal the 
spectral structures in it as good as possible. For a given 
material, the values are not necessarily the same in both 
panels. 




10 40 
Wavelength [/im] 




60 70 80 

Wavelength [/im] 



spectrum suggests that the bulk of the material around 
IRAS 16342 must be cold. A low temperature may also be 
implied by the 69 (im feature of forsterite. This intrinsi- 
cally weak feature is detected either when the forsterite 
is cold or very abundant. The feature at 4 8 ixm has been 
found in the spectra of many sources (e.g. Molster et al. 
2002a). Its nature is unclear at this time. FcSi has been 
proposed as a possible carrier l|Ferrarotti et al.1 l2000l) 
but th is identification remains debatable l|Molster et alJ 
2002a). The 48 /im feature is seen in emission, superposed 



Table 2. A list of the wavelength positions of the spec- 
tral features found in the 2 — 200 /xm spectrum of IRAS 
16342-3814. Where available an identification of a feature 
is given. Note that the amorphous silicates are positioned 
at approximately 10.9 /xm, which is at a longer wavelength 
than the normal 9.7 /xm. This difference is probably due to 
saturation of the spectrum in the 10 /jm region, making it 
difficult to determine a suitable continuum in this region 
and affecting the shape and position of the feature in the 
continuum subtracted spectrum. 



A(/xm) 


Identification 


Remarks 


10.9 ±0.4 


Amorphous Silicates 


absorption 


18.5 ±0.1 


Amorphous Silicates± 


absorption 




Forsterite±Clino-Enstatite^ (?) 


absorption 


23.5 ±0.1 


Forsterite±Clino-Enstatite t (?) 


absorption 


27.6 ±0.1 


Forsterite±Clino-Enstatite t (?) 


absorption 


29.8 ±0.1 


Diopside(?) 


absorption 


30.6 ±0.1 


? (Unidentified) 


absorption 


33.6 ±0.1 


Forsterite 


absorption 


36.1 ±0.1 


? (Unidentified) 


absorption 


36.9 ±0.1 


? (Unidentified) 


absorption 


45.3 ±0.5 


Crystalline Water Ice± 


absorption 




Crystalline Silicates(?) 




48.0 ±0.1 


? (Unidentified/ Artifact?) 


emission 


62.7 ±0.3 


Crystalline Water Ice± 


emission 




Diopside 


broad 


63.1 ±0.1 


[OI] 


emission 


68.8 ±0.1 


Forsterite 


emission 


158.0 ±0.1 


[CII] 


emission 



' By analogy with HD 161796 we tentatively identify clino- 
enstatite in the spectrum of IRAS 16342. 



on the 45 /xm feature of crystalline water ice, which is in 
absorption. If real, this suggests that the carrier of the 
feature has a different spatial distribution than the crys- 
talline water ice, otherwise both materials would be seen 
in absorption. An alternative explanation is that the fea- 
ture is an artifact in the responsivity function of the ISO 
LWS detector. More study will be needed to explore this 
issue. 



3.3. The Mass-loss history of IRAS 16342 

No other OH/IR star known to us displays a spectrum 
reddened to the extent as is seen in IRAS 16342. Also, no 
other OH/IR star is known to have crystalline silicate fea- 
tures in absorp t ion at wavelengths up to 45 /im. Though 
ISvlvester et al.l l)l999|) showed that crystalline silicate ab- 
sorption features are present in OH/IR stars, they did not 
detect them longwards of 25 /im. The behaviour of the 
spectrum and the crystalline features in it is the result of 
a very high optical depth in the dust shell, which suggests 
that the star must have an extremely large mass-loss rate. 



6 



C. Dijkstra et al.: The mineralogy, geometry and mass- loss history of IRAS 16342-3814 



We modeled the spectrum of IRA S 16342 using th e 
radiative transfer programme modust l)BouwmarJl200l() . 
We did not try to make a detailed fit of the spectrum, 
but rather wanted to get an order of magnitude estimate 
on the mass-loss rate. We assumed an M9III AGB cen- 
tral star (T e ff~2670 K and R+ = 372 R Q ) surrounded by 
a spherical circumstellar env elope (compos e d of o livine 
and iron particles, following iKemper et al ] 112002m cre- 
ated by mass-loss from the central star. We derived 
M w 10" 3 Moyr" 1 . In this case the best match between 
the peak position of the synthetic spectrum and that of 
the data was obtained. Also, the lack of near infrared 
(2 — 10 /mi) flux was best modeled in this case. 

The crystalline silicate absorption features in the ISO 
spectrum require that the envelope is optically thick up to 
about 45 (Ltm. This implies that the lack of near infrared 
flux must be due to a large opacity towards the central 
star, contrary to other post-AGB objects where the lack 
of near in frared flux is due to th e lack of dust close to the 
star (e.g. Irloogzaad et al.ll2002|) . The high mass-loss rate 
suggested by our preliminary model is consistent with this 
large opacity, and can thus explain the lack of near in- 
frared flux. However, more detailed modelling in the form 
of 2D radiative transfer calculations is needed in order to 
properly constrain the mass-loss rate of IRAS 16342. 

The high mass-loss discussed above is only expected 
to occur during the SW phase of the AGB. Combined 
with the proto-PN nature of IRAS 16342, we conclude 
that IRAS 16342 must be a very young proto-PN, per- 
haps the youngest one observed so far. The star seems to 
be right at the transition point between the AGB phase 
(which explains its extreme OH/IR type spectrum) and 
the proto-PN phase. This makes IRAS 16342 a very in- 
teresting object for the study of the transition between 
OH/IR type AGB stars and proto-PN or post-AGB stars. 
The young proto-PN nature of IRAS 16342 is in agreement 
with SAH99. Still, the possibility of a more evolved central 
star can not be excluded, since the H2O maser velocities 
in the lobes (see Sect. ^ are more in favour of an A or 
late type B star. It is difficult to understand how a later 
type star can drive such a fast wind. 

4. The infrared images 

4.1. Description of the images 

In Fig. we present the deconvolved ISAAC and TIMMI2 
images of IRAS 16342. The images were taken in the L, 
N7.9, N9.8, N2 and Q bands. The central wavelengths of 
these bands, Ao, are given in Table The lower right 
image shows the division of the N2 band image by the Q 
band image. 

IRAS 16342 reveals a bipolar morphology in the L 
band. The bipolar nature is similar to that seen in the 
optical images taken with the HST by SAH99. In both 
cases we see two lobes, surrounded by a faint halo that 
follows the contours of the lobes. The TIMMI2 images re- 
veal an ellipsoidal morphology, which may represent the 



Table 3. Size d and position angle P. A. of the bipolar 
nebula in IRAS 16342 for different bands/ wavelengths. As 
a function of wavelength, the size remains approximately 
the same, while the position angle increases. For complete- 
ness we have also listed the V band values, taken from 
SAH99. 



Band 




d(") 


P.A.( °) 


V 


0.55 


2.6 ±0.3 


71 ±2 


L 


3.78 


3.2 ±0.2 


64 ±1 


N7.9 


7.9 


1.3 ±0.3 


74 : 


N9.8 


9.8 


1.8 ±0.3 


76 ±2 


N2 


10.6 


1.4 ±0.3 


51 ±2 


Q 


20.0 


1.4 ±0.3 


51 ±2 



bipolar nebula in the ISAAC and HST images, but at 
lower spatial resolution. Some caution is needed here 
however. The optical seeing during the observa- 
tions with ISAAC and TIMMI2 were 1.2-2.0" and 
0.9-2.0" respectively (see Sect. Inland Sect. l23jl . 
In general, the seeing is better at longer wave- 
lengths under the same optical seeing. Therefore 
we might also expect the TIMMI2 images to have 
an equal or even better spatial resolution than the 
ISAAC image, in which case differences seen be- 
tween the images are intrinsic to IRAS 16342. 

For each image, the size d of the nebula along its ma- 
jor axis is listed in Table together with P.A., its posi- 
tion angle. We define the position angle as the angle be- 
tween the major axis of the bipolar nebula and the north- 
south line, where 0° represents north and 90° represents 
east. Contrary to the size of the nebula, P. A. changes 
significantly as a function of wavelength, showing a de- 
crease with increasing wavelength. For the TIMMI2 im- 
ages, both d and P. A. were measured by fitting two di- 
mensional Gaussians to the images. The Full Width at 
Half Maximum (FWHM) of the Gaussian along its ma- 
jor axis and its tilt angle (i.e. P. A. — 90°), were used to 
obtain d and P. A. respectively. For the ISAAC image we 
estimated d and P. A. by eye, while for the HST image we 
either used values given by SAH99 or derived these values 
from their data. 

The N7.9 band and the Q band reveal (sub)structure 
additional to the ellipsoidal structure mentioned in the 
previous paragraph. The N7.9 band shows some extra ex- 
tended emission in the north-east direction. Also, com- 
pared to the other bands, the morphology in the N7.9 
band seems to be less smooth, with more spatial bright- 
ness variations. However, due to the low signal-to-noise 
quality of the data it is difficult to determine if the struc- 
tures are real. Therefore, they should be treated with cau- 
tion. The Q band image reveals extended emission, or a 
halo, surrounding the central ellips. The halo is roughly 
spherically symmetric, with a gap in the north eastern di- 
rection, and approximately 6" in diameter, i.e. twice the 
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Fig. 3. The deconvolved L band (ISAAC) and N7.9, N9.8, N2 and Q band (TIMMI2) images of IRAS 16342. The 
lower right image shows the division of the N2 band image by the Q band image. Each image measures 8" x 8", with 
north is down and east is right. The contours indicate levels of equal surface brightness. No attempts were made to 
give absolute flux values (also see Sect. 12.31 andl4")l. In order to enhance the contrast of the images they were multiplied 
with a power law (of power 4). More details and a discussion of the images are given in Sect. 




cles usually scatter poorly at wavelengths large compared 
to their radii. This very red source may either be the cen- 
tral star of IRAS 16342 or the inner regions of the torus. 
If it is the central star, this suggests that IRAS 16342 only 
recently entered the post-AGB phase, in agreement with 
our results from Sect. 13.31 

An alternative to the scattering scenario in the L band 
might be that the observed light is caused by thermal dust 
emission. This is unlikely however. The dust temperature 
in the lobes at a distance r from the central star, Xd, may 
be represented by Td = T c g(r / R+) °' 4 . This expression 
assumes that the lobes are optically thin for photospheric 
radiation and that the dust emission efficiency scales as 
Q(X) oc A -1 . Using Wiens law, we need Td^650 K to 
emit efficiently at 3.78 fim. With T cff ~2670 K and R* = 
372 R Q (also see Sect. 13 .3|) we find that the radius of the 
region that can be heated to this temperature is r~ 60 AU. 
Using the distance of 2 kpc given by SAH99 and an an- 
gular diameter of 3", the observed radius of the lobes is 
3000 AU. This is 50 times larger than the expected typi- 
cal size of the emitting region derived above. We therefore 
exclude thermal dust emission as the source of radiation 
in the ISAAC image. Non equilibrium emission from small 



size of the bipolar nebula. The halo can also be seen in 
the lower right part of Fig. [21 displaying the division of 
the N2 band and the Q band images. The division shows 
that the halo is not uniformly bright. 

In ge neral our infrared im ages nicely agree with those 
made bv lMeixner et alJ l)l999|) (also see Sect.QJ. The halo 
we claim to see at 20.0 /mi has not been reported by 
iMeixner et al.l l|l999|) however. It will therefore be inter- 
esting to make additional observations at 20.0 /im in the 
future. 

4.2. The morphology of IRAS 16342 
4.2.1. The lobes 

Given the similarities with the HST image, it seems nat- 
ural to explain the bipolar nebula in the L band as due 
to reflection as well. Using this assumption, we will show 
in Sect. 14.31 that the dust particles in the lobes must be 
small (<0.1 /itm). If scattering by small particles is in- 
deed responsible for the presence of the bipolar nebula 
in the ISAAC image at 3.78 /xm, a very red bright source 
is needed to illuminate the nebula, since small dust parti- 
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grains is also unlikely, since t his requires the presence of 
UV photons l)Andriesselll978|) and these are not likely to 
be present in the system. However, it can not be ruled out 
completely. 

The behaviour of P. A. as a function of wavelength sug- 
gests that (as expected) the dominant source of radiation 
changes from scattering to thermal emission. At the short- 
est wavelengths the source is closely aligned in the direc- 
tion of the bipolar scattering nebulae seen in the HST im- 
ages; at longer wavelengths the source tends to align more 
with the torus described in Sect.^ which is expected to be 
dominated by thermal emission. Still, even up to 20 fim 
the elongation is approximately along the major axis of 
the nebula seen with HST and ISAAC, suggesting that 
the torus must be very cold (< 150 K). This is in agree- 
ment with the very red ISO spectrum, which is most likely 
dominated by thermal emission from the material in the 
torus. 

4.2.2. The halo 

We expect the halo in the Q band image to be due to 
a previous, more spherical, episode of mass-loss on the 
AGB. Given the size of the halo (6") and the bipolar neb- 
ula (3") and assuming an outflow velocity of 15 km s _1 
during the AGB phase, and a distance of 2 kpc (SAH99), 
we find that the mass in the halo was ejected between 
approximately 950 and 1900 years ago. The Q band im- 
age therefore suggests that the transition from a (more) 
spherical to a more axial symmetric form of mass-loss for 
IRAS 16342 only occured near the very end of the AGB. 
We speculate that the non uniform brightness distribution 
may indicate that there are density enhancements in the 
halo, but this remains questionable. More observations are 
required to investigate this possibility. 

4.3. Limits on grain sizes 

Using a simple model, we can put constraints on the par- 
ticle size of the dust in both the torus and the lobes. Now 
we will explain the model in detail and discuss its results. 
We will show that there are two independent methods that 
allow us to derive the optical depth in the torus along the 
line-of-sight towards the eastern lobe. The first method 
depends on the dust particle size in the lobes, while the 
second method depends on the particle size in the torus. 
Both methods must yield the same optical depth for a 
given wavelength. Applying this constraint and calculat- 
ing the optical depth with both methods for three different 
wavelengths, yields the particle size in both the lobes and 
the torus. 

Our model makes use of the following assumptions. 
The eastern lobe in the HST and ISAAC images is signif- 
icantly redder than the western lobe. In order to explain 
this additional reddening, we assume the eastern lobe to 
be located completely behind the torus, and the western 
lobe to be positioned completely in front of the torus (see 



Fig. 4. The assumed morphology for IRAS 16342 used to 
determine the dust particle sizes in both the torus and the 
lobes. For details and a discussion see Sect. 14.31 




t 

/ e ,A = /oGA(180-i,ai)e-^ 

\f W Direction to Earth 

Iw,x = IoG\(i, ai) 

Fig. 0} . The torus acts as a screen that causes additional 
reddening of the eastern lobe with respect to the western 
lobe. We assume that scattering is the dominant source 
of emission in the lobes, and that both scattering lobes 
are intrinsically equally bright. Even in the absence of the 
(extinction due to the) torus, the lobes will not appear 
equally bright as seen from Earth, due to the orientation 
of the system. The orientation causes the eastern lobe to 
be dominated by back scattering, while the western lobe 
is dominated by forward scattering. As a result, the east- 
ern lobe will appear less bright than the western lobe. 
Interstellar reddening is assumed to be the same for both 
lobes, and scattering of stellar radiation within the torus is 
assumed to be unimportant (due to its large radial optical 
depth). 

In the first method to derive the optical depth in the 
torus we measure the intensities of the eastern and western 
lobes, I e ^\ and I w ,\, as seen from Earth (see Fig.^J. These 
are given by 

h,x =/ GA(180-z, ai )e~ rA (1) 

I w ,\=hG x {i,ai) (2) 

respectively, where Iq is the intrinsic intensity of each 
lobe. a\ is the characteristic particle size in the lobes, 
while G\(9,a\) is the scattering phase function of a par- 
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ticle with radius a\ irradiated at wavelength A and scat- 
tering at an angle 9. The phase function was calculated 
for size distributions of sp herical grains using Mie theory 
l|Bohren fc Huffmanlll98^ . Finally, i is the inclination an- 
gle of the system, defined as the angle between the major 
axis of the nebula and the line-of-sight. We take i = 40° 
(SAH99). t\ is the optical depth of the torus along the 
line-of-sight towards the eastern lobe. From Eq. JIJ and 
J5J) we find 



t\ = In 



I w ,xGx(lS0-i,ai) 
Ie,\G\(i,a\) 



(3) 



This yields t\ for different values of oi. 

The intensity ratio I w .\/Ie.\ was determined from the 
magnitude difference between the lobes, Am, which was 
either given by SAH99 or we measured it ourselves. This 
implies that we replaced the intensity ratios by flux ratios. 
The values for the magnitude differences are listed in Table 
01 together with the corresponding intensity ratios. 

We describe the grain size distribution used in the 
calculation of G\(9,a\) by a power law rt(a)oca _m with 
m = —3. This value was chosen for calculation conve- 
niences. Differences with mode ls using the more g enerally 
adopted ISM value of -3.5 fe.g. lMathis et alll977|) are ex- 
pected to be small. The effective width of the size distribu- 
tion was kept fixed. The grains in the lobes were assumed 
to be composed of silicates, for which we used the optical 
constants of 'cold' silicates derived bv ISuhl l)l999|) . Using 
the op tical constants of 'warm' silicates derived bv ISuhl 
(1999) yields very similar results (not shown). 

In the second method we derive t\ from the well known 
relation 



t\ = n d 7ratQ^ xt (a t )£ 



(4) 



where rid is the dust particle number density in a dust col- 
umn of length L along the line-of-sight; at and Q^ xt (at) 
are the characteristic radius and extinction efficiency of 
the particles in the torus. For the calculation of the ex- 
tinction efficiency we adopt the same dust properties as 
those assumed above for the calculation of the scattering 
phase functions. 

Eq. J3J and @ should yield the same result for t\ 
at a given wavelength. Ratios between optical depths are 
independent of rid and L and therefore allow a relatively 
easy derivation of the particle sizes in the lobes and the 
torus. 

Using Eq. J2J| and the appropriate intensity ratios (see 
Table 0} and phase functions, we measured t\ at 0.55, 
0.79 and 3.78 /im for different assumptions for the char- 
acteristic particle size in the lobes, a\. In Fig. is 
plotted versus with a\ as a free parameter (see line 
indicated with 'I-ratios'). We also calculated Q x xt (a t ) at 
the above wavelengths for different values of the charac- 
teristic particle size in the torus, at. Since T\ccQ x xt (at) , 
we must have 



TO. 55 _ Qo X 55(°t) 



r 0.80 



?0 X 8o( fl t) 



(5) 



Fig. 5. IiLM - against Is ^ m - for different particle sizes in 

° To SO ° T3.78 * 

the torus (labeled with 'Q-ratios') and different particle 
sizes in the lobes (labeled with 'I-ratios') of IRAS 16342. 
Also shown are the p ositions for the ISM (based on 
ISavaee fc Mathislll979|) and a grey absorber (indicated by 
stars) . The arrows indicate directions of increasing particle 
sizes. For some points (indicated by diamonds) the max- 
imum particle size of the used size distribution is shown 
(in /mi) . For a detailed explanation and discussion of this 
figure see Sect. 14.31 



w 



Grey Absorber ,-~ 2 
10.0 




Direction of increasing particle radius.a 

10.0 '00.0 



Table 4. The intensity ratios, I w ,\/Ie.\, at given wave- 
lengths. This quantity is relevant for the calculation of t\ 
in Eq. ©. Listed are the wavelength; the magnitude differ- 
ence between the lobes; the intensity ratios. The adopted 
errors on the magnitudes are expected upper limits. 



Ao(amti) 


Am 


Iu,,X 
Ie.X 


0.55 


2.8 ±0.1 


13.2±0.5 


0.79 


2.2 ± 0.1 


7.6±0.3 


3.78 


0.64 ±0.1 


1.80±0.07 



and 



T 0.80 



_ Qo X 8o( a t) 

T 3.78 Qf^si ^) 



(6) 



independent of n d and L. In Fig. we again plot 
versus ^^-j this time with at as a free parameter (see 
line indicated with 'Q-ratios'). The intersection of the two 
curves provides a consistent set for a\ and at. 

For the lobes we find that we have a size distribution 
with a maximum particle size aj nax = 0.09±0.02 /im, while 
for the torus we find a{° ax = 1.3±0.2 /im. Both values 
are independent of the adopted effective width of the size 
distribution, which we take to be the same for particles 
in both the torus and the lobes. A test with amorphous 
olivines instead of the cold 'Suh' silicates yielded a similar 
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result. The errors on the particle sizes are derived from 
an error analysis on Eq. J2J taking the errors on Am and 
I w ,\/Ie,\ into account for each wavelength (see Table^J. 
We are interested in the maximum particle size of the 
distribution instead of the typical (i.e. effective) particle 
size. This is because we want to set upper limits on the 
particle sizes found in the lobes and torus. The typical size 
simply scales linearly with the maximum particle size. 

We thus derive that the maximum particle size in the 
torus is about an order of magnitude larger than the max- 
imum particle size in the lobes. In an AGB outflow, the 
particle size of the grains formed in the outflow is pro- 
portional to the density of the e nvelope at the dust con- 
densation radius l|Dominiklll992|) . Larger densities provide 
larger grains. The spatial distribution of large particles in 
a torus geometry suggests that this material is expelled 
from the AGB star in an aspherical high density wind. 
In a natural way, this accounts for the presence of large 
grains in the torus and small grains in the lobes. However, 
the bipolar morphology of IRAS 16342 may also have orig- 
inated from a later phase, i.e. after the AGB. In this case, 
the AGB wind may have been spherical, but due to the 
interaction of this wind with the (more recent) fast bipolar 
molecular outflow, the morphology may have changed into 
the observed torus like spatial distribution. The interac- 
tion of these two flows may have caused a shattering of the 
dust particles in the lobes into smaller ones, while leaving 
the grains in the torus intact. Regardless of the correct sce- 
nario, the presence of the bipolar nebula in IRAS 16342 at 
optical wavelengths suggests that the mass-loss rate near 
the central star very recently must have dropped, other- 
wise the optical photons in this nebula would not be able 
to escape. Obviously the star has already managed to cre- 
ate a bipolar cavity while the optical depth of the dust 
shell in the equatorial torus is still huge. 

The derived particle size for the torus is large com- 
pared to the typical ISM value s, which are of the order 
of a ~ 0.01 to ~ 0.2/zm (e.g. iMathis et aflll977|) . but 
in agreement with other objects for which large mass- 
loss rates are expected. Exam ples of such source s are 
the binary system AFGL 410 6 llMolster et aljll999h and 
the post-AGB star HD161796 i|Hoogzaad et all2002|) . The 
large grain size in the torus explains why the dust parti- 
cles in the torus of IRAS 16342 are grey compared to the 
ISM. 

We finally note that when the emission in the ISAAC 
image at 3.78 /im is caused by non equilibrium emission 
from small dust grains, as was considered in Sect. 14.2.11 
the resulting particle sizes in the lobes and torus are very 
similar to those derived above. The possibility of non equi- 
librium emission by small grains at 3.78 /zm can simply be 
modeled by leaving ^3.78(180 — i, a{) and i r 3.7s(i, ai) out of 
equation|31for T3.78, i.e. by assuming that at 3.78 /zm there 
is no scattering. 

5. Conclusions 



1. The ISO spectrum of IRAS 16342 is that of an ex- 
tremely reddened OH/IR star. In addition to amor- 
phous silicate absorption features at 10 and 20 /im, 
the spectrum shows crystalline silicate features in ab- 
sorption at wavelengths up to 45 /im. Only at longer 
wavelengths the crystalline features are in emission. No 
other OH/IR star is known to show this behaviour. We 
estimate the mass-loss rate of IRAS 16342 may be as 
large as 10~ 3 Moyr^ 1 . This high rate explains the lack 
of flux in the near-IR as well as the peak wavelength 
in the SED. 

2. We identify forsterite (Mg 2 Si04), diopside 
(MgCaSi20e), and possibly clino-enstatite (MgSiOs) 
in the spectrum. Besides these crystalline silicates, 
also crystalline water ice is detected. All these mate- 
rials are also found around other (high M), oxygen 
rich, evolved stars. At 48 /Ltm a (to date) unidentified 
feature is seen. 

3. IRAS 16342 must be a very young proto-PN, perhaps 
the youngest one observed until present. The central 
star seems to be right at the transition point between 
the AGB phase and the proto-PN phase. 

4. The deconvolved L band ISAAC image of IRAS 16342 
at 3.78 fim reveals a bipolar scattering nebula. Both 
lobes of the nebula are separated by a dark equato- 
rial waist, or torus. This morphology is consistent with 
that found in HST images. The presence of the scat- 
tering nebula in the L band suggests that the central 
star of IRAS 16342 may be very red, which would be 
in agreement with the young proto-PN nature of the 
system. The TIMMI2 Q band image reveals thermal 
emission from an extended, roughly spherical halo of 
~ 6" in diameter. We expect the halo to be due to 
a previous, more spherical, phase of mass- loss on the 
AGB. This halo material was ejected between approx- 
imately 950 and 1900 years ago. 

5. TIMMI2 images at wavelengths intermediate to the 
L and Q band reveal an ellipsoidal morphology, which 
may represent the bipolar nebula but at lower spatial 
resolution. The position angle of the nebula changes as 
a function of wavelength, suggesting that the nature 
of the radiation changes from scattered light to thermal 
emission (by cold material). Still, even up to 20 /im the 
elongation is approximately along the major axis of 
the nebula seen in the optical and near-IR, suggesting 
that the equatorial torus must be very cold. This is in 
agreement with the very red ISO spectrum, which is 
most likely dominated by thermal emission from this 
torus. 

6. Using a simple model for the geometry of IRAS 16342, 
we calculate the dust particle sizes in both the torus 
and the lobes. The grains in the lobes have a maximum 
particle size ap ax = 0.09 /im, while for the torus we 
find a™ ax = 1.3 /im. The relatively large torus grains 
suggest that the bipolar morphology of IRAS 16342 
may already have been created during the AGB phase. 



The results of our study can be summarised as follows: 
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